Crossing intervals of non-Markovian Gaussian processes 
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We review the properties of time intervals between the crossings at a level M of a smooth station- 
ary Gaussian temporal signal. The distribution of these intervals and the persistence are derived 
within the Independent Interval Approximation (II A). These results grant access to the distribution 
of extrema of a general Gaussian process. Exact results are obtained for the persistence exponents 
and the crossing interval distributions, in the limit of large \M\. In addition, the small time behav- 
ior of the interval distributions and the persistence is calculated analytically, for any M. The IIA 
is found to reproduce most of these exact results and its accuracy is also illustrated by extensive 
numerical simulations applied to non-Markovian Gaussian processes appearing in various physical 
contexts. 

PACS numbers: 05.40.Fb, 02.50.Cw, 02.50.Ey 



I. INTRODUCTION 

The persistence of a temporal signal X is the probabil- 
ity that X{t') remains below (or above) a given level M, 



for all times t' £ [0, i\. This problem has elicited a large 



body of work by mathematicians 
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physicists, both theorists 10 
|l|, P 1^, |l| and 

Persistence properties have been measured in as different 
systems as breath figures p^ , liquid crystals laser- 
polarized Xe gas p4| , fluctuating steps on a Si surface 
p5[ , or soap bubbles p6[ . 

Although the persistence is a very natural and easy 
quantity to define ~ in a sense, its physical meaning is 
easier to explain to the layman than the meaning of a 
two-point correlation function - this quantity turns out, 
in practice, to be extremely complex to deal with ana- 
lytically. In fact, exact results have been obtained in a 
very limited number of cases, as far as non-Markovian 
processes are concerned |§, ^, ■ 

The mathematical literature has mainly focused on 
evaluating the persistence 



P<{t) = Vioh{X{t') < M, t' e [0, t]), 



(1) 



mostly for Gaussian processes {i.e. processes X{t) for 
which the joint distribution of X{ti), ...,X{tn) is Gaus- 
sian) and for large \M\, a regime where efficient bounds 
or equivalent have been obtained ||, ||] . Recently (D , 
and for Gaussian processes only, a numerical method to 
obtain valuable numerical bounds has been extended to 
all values of M, although the required numerical effort 
can become quite considerable for large t. 

Physicists have also concentrated their attention to 
Gaussian processes |l^, ^ ^ |l5|, |l^, , which are often 
a good or exact description of actual physical processes. 
For instance, the total magnetization in a spin system 



]l6| , [r^ , the density field of diffusing particles |15|, or the 
height profile of certain fluctuating interfaces |19, |4| 
are true temporal Gaussian processes. Two general meth- 
ods have been developed, focusing on the case M = 0, 
which applies to many physical situations. The first one 
|]l2| , [l3| , [l^ is a perturbation of the considered process 
around the Markovian Gaussian process, which has been 
extended to small values of |A/| |1^. Within this method, 
only the large time asymptotics of P< (t) is known, lead- 
ing to the definition of the persistence exponent (see be- 
low). The alternative method, using the Independent In- 
terval Approximation (IIA) |lj, gives very accurate 
results for "smooth" processes, that is, processes having a 
continuous velocity. Initially, the IIA remained restricted 
to the case M = |lq| , but it has been recently general- 
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ized to an arbitrary level M |20 . 

The persistence probability is also intimately related 
to another important physical quantity: the probability 
distribution of the extrema of the considered process X . 
For instance, the quantity P<(t) can be also viewed as 
the probability that the maximum of X in the interval 
[0, t] has remained below the level M. Thus, the distribu- 
tion of the maximum of X is simply the derivative with 
respect to M of P<{t). The distribution of the extrema 
has been analytically obtained for the Brownian process 
p7| , p8[ , but its derivation remains a formidable task for 
general non-Markovian Gaussian processes. On this ac- 
count, the persistence problem is also related to extreme 
value statistics which is a notoriously difficult 

problem for correlated variables, and which has recently 
attracted a lot of attention among physicists . 

Hence, the persistence problem has obvious applica- 
tions in many other applied and experimental sciences, 
where one has to deal with data analysis of complex sta- 
tistical signals. For instance, statistical bounds of noisy 
signals are extremely useful for image processing (for in- 
stance in medical imaging or astrophysics |^2| ), in order 
to obtain cleaner images by correcting spurious bright or 
dark pixels (l| ^ . In general, it is important to be able to 
evaluate the maximum of a correlated temporal or spa- 
tial signal originating from experimental noise. The same 
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question can arise when the signal lives in a more abstract 
space. For instance, in the context of genetic cartogra- 
phy, statistical methods to evaluate the maximum of a 
complex signal has been exploited to identify putative 
quantitative trait loci Jsst - Finally, this same problem 
arises in econophysics or finance, where the probability 
for a generally strongly correlated financial signal to re- 
main below or above a certain level is clearly of great 
interest Q. 

In the present paper, we are interested in the persis- 
tence of a non-Markovian Gaussian process, which can 
be either stationary or scale- invariant. More specifically, 
this study focuses on the properties of the distributions 
of time intervals during which the considered process re- 
mains below or above a given level M . We shall see that 
these distributions are simply related to the persistence 
itself, and contain valuable physical information. 

We now summarize the content of the present work, 
and mention briefly our main results. In section ||, we in- 
troduce our main quantities of interest, and among them: 
«. the persistence P>{t) and P<{t) which measure the 
probability to remain above or below the level M up to 
time t] ii. the associated distributions of ± time intervals, 
P± (t) , during which the process remains above or below 
the level M. For a general level M, we then shortly re- 
view the IIA calculation introduced in , which leads to 
analytical approximate expressions for the Laplace trans- 
form of P±{t), P>{t), and P<{t). In fact, we present a 
simpler formulation compared to ||2^ , as well as new ad- 
ditional results, which permit a fast and efficient numer- 
ical implementation of the IIA results. In section HI, 
we first introduce several examples of physically rele- 
vant smooth Gaussian processes, on which the general 
results of this paper are numerically tested. Within the 
IIA, we then present the calculation of the persistence 
exponents 0±{M), defined, for stationary processes, by 
P±{t) ^ ex-p{—6±t), when t +00. We also obtain 
exact estimates for 9±{M), for large M > 0, including 
exact bounds for P>{t) and 9+{M), and compute the ex- 
act asymptotic distributions P±{t). For large M, P+{t) 
takes the form of the Wigner distribution, whereas P_ (t) 
becomes Poissonian. All these results are correctly re- 
produced by the IIA, except for the large M asymptotics 
of 0+{M). Finally, in section [V, we obtain the exact 
small time behavior of P± (t) , P> (t) , and P< (i) , for "very 
smooth" processes (a term to be defined in the next sec- 
tion H), and find that the IIA again reproduces these 
exact results. For marginally smooth processes, we also 
obtain exact results, but the IIA is not exact in this case, 
although its results remain qualitatively correct and even 
quantitatively accurate. Many of these results are also 
obtained by means of simple heuristic arguments, help- 
ing to understand them on physical grounds. This study 
is also illustrated by means of extensive numerical sim- 
ulations, revealing a very satisfying accuracy of the IIA 
for moderate M, in a regime where, contrary to the limit 
M ±00, there are no available exact results for P±{t), 
P>{t), and P<{t). 



II. INDEPENDENT INTERVAL 
APPROXIMATION (IIA) 

In this section, we introduce our main physical quanti- 
ties of interest - interval distributions, persistence, sign 
autocorrelation, constrained densities of crossings - and 
relate their general properties. We then summarize the 
IIA calculation introduced in [20|] , and obtain more ex- 
plicit expressions for the interval distributions and the 
persistence probability, in the case of a Gaussian process. 
These new results will prove useful in the next sections 
more specifically devoted to the interval distributions. 



A. Introductory material and notations 

One considers a stationary non-Markovian Gaussian 
process X{t) of zero mean and unit variance. Its distri- 
bution at any time is then 



gix) 



e 2 



2n 



and we define its cumulative distribution as, 
G{X)= j g{x)dx = l-G{X). 

J —00 



(2) 



(3) 



Due to its Gaussian nature, such a stationary process 
is entirely characterized by its two-point correlation func- 
tion. 



(4) 



It is understood that the process starts from t = — 00, 
so that the distribution of the initial condition X{0) at 
t = is given by Eq. (||). In addition, all derivatives of 
X{0), when they exists, are random Gaussian variables 
of zero average and second moment 



X("'(0)1 ) = (-l)"/(2n)(0), 



(5) 



where the superscript (n) refers to a derivative of order 
n. 

The process is assumed to be "smooth" , although this 
constraint will be sometimes relaxed in the present paper. 
By "smooth" , we mean that the velocity of the process is 
a continuous function of time. "Very smooth" processes 
will have a differentiable velocity. In particular, the most 
general stationary Markovian Gaussian process defined 
by the equation of motion. 



dX 
~dt 



= -\X + ^/2\^ 



(6) 



does not belong to this class of smooth processes {ri{t) is a 
Gaussian 5-correlated white noise). In practice, a process 
is smooth if its correlator f{t) is twice differentiable at 
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t = 0, and is very smooth if the fourth derivative of f{t) 
exists at t = 0. The process introduced in Eq. has a 
correlator 



/(0=exp(-A|i|), 



(7) 



which has a cusp at t = 0, and is thus not even differen- 
tiable at t = 0. 

Throughout this paper, we will be interested in the 
probability that the process X remains below or above 
a certain threshold M, for all times in the interval [0,t]. 
In particular, the smoothness of the process ensures that 
the M-crossings {i.e. the times for which X{t) = M) are 
well separated, with a finite mean separation between 
them, denoted by r. For the Markovian process men- 
tioned above, the Af-crossings are distributed on a frac- 
tal set of dimension 1/2, and the mean interval is then 
T = 0. For smooth process, r can be computed by eval- 
uating the mean number of M-crossings during a time 
interval of length t, 



N{t) 



\X'{t')\5{X{t')- M)dt' 



(8) 



The correlation functions between the position X and the 
velocity X' is 



{x{t)x\t')) = -nt-t'). 



(9) 



By time-reversal symmetry, and since f{t) is twice differ- 
entiable at t = 0, one has 



{X{t)X\t))^-nQ) = 0, 



(10) 



so that the position and the velocity are uncorrelated at 
equal time. The distribution of the velocity X'{t) is a 
Gaussian of zero mean and second moment 



so that 



(X'2(t)> =-/"(0)=a2. 



\X'{t)\) = f-^. 



Finally, N{t) = {\X' {t)\)g{M)y.t, which leads to 



(11) 



(12) 



(13) 



We also define the distributions of time intervals be- 
tween M-crossings, P+{t) and P-{t), during which the 
process remains respectively above and below the level 
M. The means of these two kinds of intervals are defined 
by T± = J^°° t P± (t) dt, and are related to r. Indeed, 
since there are as many -I- and — intervals, as they sim- 
ply alternate, one has 



T = — = e 2 



In addition, r_|_ / t_ is also equal to the ratio of the times 
spent by the process X above and below the level M, 
i.e., G{M)/G{M). Finally, we obtain 



2tG(M), 



2tG(M). 



(15) 



We now introduce the persistence of the process X, de- 
fined as the probability that it does not cross the level M 
during the time interval [0,t]. More precisely, we define 
P>{t) and P<{t) as the persistence, knowing that the 
process started respectively above and below the level 
M. In other word, P>{t) is also the probability that the 
process remains above the level M during the considered 
time interval, and P<{t) is the probability that X re- 
mains below the threshold M. Note that the persistence 
probes the entire history of the process, and is therefore a 
infinite-points correlation function: for instance P< (t) is 
the probability that the process remains below the level 
M between the times and dt, and dt and 2dt,... and 
t — dt and t. For non-Markovian processes, it is thus un- 
derstandable that this quantity is extremely difficult to 
treat analytically, and there are very few examples where 
P> (t) or P< (t) can be actually computed H, |^, ||] . 

P> (t) and P< {t) are intimately related to the interval 
distributions P+(<) and P-{t) by the relations j2^, 



/ + OC 
{t'-t)p+{t') 



dt' 



(16) 



P<(t) = tZ^ j {t' -t)P_{t')dt'. (17) 

Indeed, if X{t) has remained below the level M up to 
time t, it belongs to a — interval of duration t' > t, start- 
ing at an initial position uniformly distributed between 
and t' — t, an interpretation which leads to Eq. (p7|). The 
above expressions can also be differentiated twice, giving 



P+{t) = T+P^{t), 

P^{t) = r^P'^it). 



(18) 
(19) 



If M = by symmetry of the process under the 
transformation X —X, one has P>{t) — P<{t) and 
P+{t) = P-(t). In general, we have the symmetrical 
relations, 

P>(i, M) = P<(t, -M), P+{t, M) = P_{t, -M), (20) 

so that we will restrain ourselves to the case M > 0. 

The knowledge of P> [t) and P< [t) also provides valu- 
able information on the distribution of the extrema of the 
considered process. The statistical properties of extremal 
events_2S an_ active field of research among mathemati- 

By definition, P< [t) is 



cians 



and physicists |31 



also the probability that the maximum of the process in 
a time interval of duration t is less than M, and P>{t) 
is the probability that the minimum of the process re- 
mains bigger than M for all times in [0,t]. Hence, defin- 

(14) 

ing Pmax(Af, t) and Pinin(M, i) as the distribution of the 
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maximum and minimum of the considered process, one 
has, 



(21) 

(22) 



Note finally that Eqs. (^jl7||lj,|ig|,y,g2|) are in fact valid 
for any stationary process, not necessarily Gaussian. 

Before presenting an approximate method leading to 
analytical expressions for the different quantities intro- 
duced above, we need to define two quantities which will 
prove useful in the following. We start with the autocor- 
relation function of d[M — X (t)] {6 is Heaviside's function: 
e{x) = 1, if X > 0; e{x) = 0, if a; < 0; 61(0) = 1/2), 

A>it) = (0[X{t)-Ad]e[X{Q)-M]), (23) 
A^{t) = {e[M-X{t)]e[M-X{0)]), 

= 2G{M)-1 + A>{t), (24) 

where the last relation is obtained by using 9{x) = 1 — 
9{—x). In addition, since the process is invariant under 
the transformation X ^ —X, one also has, 



^>(Af,t) = A<(-A/,<). 



(25) 



For a Gaussian process, these quantities can be explicitly 
expressed in terms of the correlation function f{t) |p5|. 



A>{t) = 
A<{t) = 



+ 00 

M 
M 



g{x)G 



g{x)G 



( M- 


- =^f{t) 




-Pit) 


' M - 






-Pit)) 



dx. 



dx. 



(26) 
(27) 



For M — 0, these integrals can be explicitly performed 
PI, giving 



1 



A> {t) ^A^{t)^- + — arcsin(/(i)). (28) 
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Finally, the time derivative of A^{t) and A^{t) can be 
simply written as 



A'^{t) = A'S) 



1 



27r 



exp 



A/2 



1 + 



(29) 

We now introduce Ny{t) and N<^{t), the mean number 
of Af-crossings in the time interval [0,t], knowing that 
the process started from X{Q) > M and X{0) < M, 
respectively. These quantities satisfy the sum rule, 

GiM) N<{t) + GiM) N>{t) = N{t) = -, (30) 

T 

which expresses the fact that a Af-crossing is either 
crossed from above, or from below M. Again, it is clear 
that 



Ny{M,t) = N<{-M,t). 



(31) 



In addition, for large time t, 

N^{t)^N^{t)^N{t) 



t +O0, 



(32) 



since the initial position of X{t) becomes irrelevant when 
t +00. On the other hand, for short time, 

7V>(i)^— , N<{t)^—, t^O, (33) 

which expresses the fact that the probability per unit 
time to meet the first A/-crossing is respectively r^^ 
and tZ^ , for + and — intervals. Note that both asymp- 
totics of Eqs. ( |3^ , ^ ) are consistent with the sum rule of 
Eq. (|30|). For a Gaussian process, N^{t) and -/V<(t) can 
be calculated after introducing the correlation matrix of 
the Gaussian vector {X{t), X{0), X'{t)), which reads 

/ 1 fit) \ 

m = fit) 1 fit) . (34) 
V fit) -no) J 

For instance, in the same spirit as Eq. (||), one has 

N^(t)=G-\M) l\\X'{t')\)^dt', 



(35) 



where (|X'(t)|)< is the average of the velocity modulus, 
knowing that X{t) ~ M, and averaged over X{0) < M: 



{\x'it)\). 



M 



dxQ 



dv 



me 2 



(27^)3/2^/ditC 
where U = {M,xq,v). Ny{t) is similarly defined as 



(36) 



N^it)=G-\M) / {\X'it')\)ydt', (37) 
Jo 



with 



{\X\t)\): 



+ 00 

dxQ 

M J -oo 



, |i;|e-iutc-iu 

dv 



(27r)3/2Vd^ 



= r-i-(|X'(t)|)<. 
|X'(t)|)< can be written more explicitly, 



(38) 
(39) 



T{\X'{t)\)<=Gib) + a 
with 

a{t) 
bit) 



--G[ah) 



e" — -T+T, (40) 



\f'it)\ 



V^2il-Pit)) 

1 - fit) 



M 



V^-Pit)-rHt)/a2' 



(41) 
(42) 



where 02 = — /"(O). Using Eq. (P9|), one finds a similar 
expression for {\X'{t)\)y, 



■(|X'(t)|)>=G(6)-a 



l-Giab) 



e" — -T+T. (43) 
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When t +00, a{t) -> and b{t) M, so that ||3|| 

(44) 



Usin g E qs. ( p5| , p7| ), one recovers the asymptotics of 
Eq. (p2|). On the other hand, when t ^ 0, we have 
a{t) ^ 1 and f{t) 1, which leads to 



(|X'W|)< ^ (|X'(i)|>> 
and we recover Eq. (ISSl). 



1 

27' 



(45) 



B. Derivation of the II A distributions 

In the previous subsection, we have introduced the 
so far unknown interval distributions P+ (t) and P_ (t) 
which are intimately related to the persistence prob- 
abilities P>(t) and F<(t) (through Eqs. (|l|,|l^,|l|,^) 
and the distribution of extrema of the process (see 
Eqs. ( |2l| , p^ )). On the other hand, for a Gaussian process, 
we have computed explicitly the autocorrelation A> (t) 
(and A^{t)) and the constrained number of M-crossings 
Ny{t) (and N^{t)) as a function of the correlator f{t) 
(see Eqs. (|29[^,||)). We will now try to relate the two 
unknown interval distributions to these two known quan- 
tities. 

Let us define P< {N, t) and P> {N, t) , as the probability 
that there are exactly N M-crossings in the interval [0,t], 
starting respectively from X{0) < M and X{0) > M. By 
definition of Ay{t), one has 

-t-00 

A>(t) = G(M)^P>(2n,t), (46) 

+00 

A<(t) = G(Af)^P<(2n,t), (47) 



since X{t) is on the same side of M as X(Q) if and only 
if the number of M-crossings in the interval [0, t] is even. 
{t) and {t) can also be simply written as 



4-00 



^>(*) = ^nP>(n,t), 

+00 

iV<(t) = ^nP<(n,t). 



(48) 
(49) 



ri=0 



Note that, by definition, one has P>(0,i) — P>{t) and 
P<(0,t) = P<(t). 

Our central approximation now consists in assum- 
ing that the interval lengths between M-crossings are 
uncorrelated jT5[ This will lead to closed rela- 

tion between (P<(iV,i), Py{N,t)) and P±{t). Using 
Eqs. (^,^,^,E9|), we will then obtain an explicit ex- 
pression of P±(t) sls a function of A>(t) (or A<(t)) and 
iV>(0 (or7V<(i)). 



Let us consider an odd value of = 2n — 1 {n > 1). 
Using the IIA, we obtain 

P<(2n-l,i) =rri/o^*iQ-(^i)x 
J*^ dt2 P+{t2 - ti) //^ dt3 P_{t^-t2)--- 

■S't2n-3 '^^2n-2 P+(^2n-2 ~ i2»i-3)x 
/i _ dt2n-l P-{t2n-l — t2,i-2)Q+{t — t2n-l), (50) 



where 



-it) 



+CX2 



P±{t')dt', 



(51) 



is the probability that a ± interval is larger than t. 
Eqs. (|5^) expresses the fact that to find 2n — 1 cross- 
ings between and t starting from X{Q) < M, we should 
find a first crossing at ti (and hence an initial — interval 
of length bigger than ti, with probability TZ^Q-{ti)), 
followed by a -I- interval of length t2 — ti, and so on, up 
to a last crossing time t2n-i, associated to a — interval 
of length t2n-i — i2n-2- Finally, there should not be any 
further crossing between t2n~i and i, hence the last fac- 
tor Q+{t — i2n-i)- AH these events have been treated 
as independent, so that there probabilities simply factor: 
this is the core assumption of the IIA. For even N — 2n 
{n > 1), one obtains a similar expression, 

P<{2n,t) = tZ^ dtiQ^{ti)x 

dt2 P+{t2 - tl) dt3 P_(t3 " ^2) • • • 



Ju 



dt 



2n-l 



P-(^2n-l — i2n-2)x 



Ito _i '^^2ri P+(^2n — t2n-l)Q-{t — t2n 



(52) 



For a given function of time F{t), one defines its 
Laplace transform, F{s) = J^°° F{t) e^''* dt. The con- 
volution products in Eqs. ( [50|j5^ ) take a much simpler 
form in the Laplace variable s: 



P<(2n-l,s) 
P<(2n,s) 



tZ'Q+Q4P+P^] 
tZ'Q^_P+[P+P-] 



n-l 



n-1 



where 



1-P±(s) 



(53) 
(54) 



(55) 



is the Laplace transform of Eq. 
conservation of probability, 



)T|). If we express the 



+ 00 



P<{t) + J2P<{N,t) = l, 



(56) 



N=l 



and after summing simple geometric series, we obtain 

L 1 - P_ (s) 



P<i-s) 



(57) 
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This relation is nothing but the Laplace transform of 
Eq. (|l7|). It is certainly satisfying, and also reassur- 
ing, that the IIA reproduces this exact relation, as well 
as the equivalent relation between P>{t) and P+{t) of 
Eq. Of course, F>(2n- lA) and P>(2n,i) satisfy 

similar equations as Eqs. ( ^ , ^2|j5^j5^ , |56| , |57| ) , obtained by 
exchanging the indices + ^ — and <<-!■>. 

Using Eqs. ( p6| , p7| , p8 49), we can now write explicitly 
the Laplace transform of the known quantities A< (t) and 
iV<(t) in terms of the Laplace transform of P±{t): 



N<is) = 



A<(s) = G(Af) 



(1 + P+)(1-P-) 



1"P+ 



N<is) 



(58) 
(59) 



Again, Ay {t) and iV> (t) are given by similar expressions 
after the substitution + — and < >, and G{M) ^ 
G(M). 

Using P4(0) = - J^°° t P±{t) dt = -T± and Eq. ^ 
one obtains the following estimates when s — > 0, 



N<{s) J, i<(s) 



G2(M) 



(60) 



The first expression in Eq. ( pO| ) is equivalent to the gen- 
eral result of Eq. (|3^), whereas the second relation ex- 
presses that for large t, A^{t) ~ G'^{M). For large s, 



N<is) 



1 



G{M) 



(61) 



which corresponds to the small time behavior of Eq. (|33| ) 
for N^{t), whereas the second relation is equivalent to 
A<{0) = G{M). 
Finally, writing 



G(M)-si<(s) 

G{M)sN<{s) 



(62) 



and using Eqs. ( |58| , |59| ), the interval distributions are 
given by 



P.{s) 



2-T_s2jV<(s)(l + f<(s)) 
2-r_s2iV<(s)(l-F<(s))' 



(63) 
(64) 



the process under the transformation X ~X, the fol- 
lowing symmetric relations hold: 

P±{-M,t) = P^{M,t), P:^{-M,t) = P<iM,t). (65) 



Using Eq. (|29|) , we find that the dimensionless function 
W{t), 



W{t) 



-rA'y{t) = -rA'^{t), 
r fit) 



a _m1 1-/ 

_ p 2 ■ T+T 

2 



exp 



1 + fit) 



, (66) 



(67) 



has a simpler analytical form than Ay(t) or A^{t). Sim- 
ilarly, we define the dimensionless auxiliary function, 
V{t), by 

rGiM)N^{t) = r{\X'{t)\)y=GiM) + V{t), (68) 
TG(M)iV^(t) = T{\X'{t)\)^=G{M)-V{t). (69) 



Using Eq. (|40|), we find that V{t) takes an explicit form 
contrary to iV> (t) or A''< (t) : 



V{t) = G{M) - G(6) -f a 



G(a 6) 



1 



Aj2 1-/ , 

e" — -1+7, (70) 



where a(t) and b(t) are simple functionals of the corre- 
lator f{t), which have been defined in Eqs. (^ 42). We 
give below the behavior of W{t) and V{t), in the limit 
t 



V{t) 
W{t) 



M 



: e 2 



■ e 2 



■fit). 



(71) 
(72) 



while one has y(0) = G(M) - 1/2 and PF(0) = 1/2. 

The Laplace transform of W(€) and V(t) can be ex- 
plicitly written 

W^(s) = G(A/) - si<(s) = G(Af) - si>(s), (73) 

and 

rG(A/)s7V>(s) = ^i^ + -i/(s), (74) 
s 

rG(M)s7V<(s) = ^li^_l/(s). (75) 



Inserting these expressions of P±{s) in Eq. (|57|), one ob- 
tains P<{s) (and P>(s)), from the sole knowledge of 
A<(f) and iV<(i) (or their Laplace transform), which 
are known explicitly for a Gaussian process. Alterna- 
tive expressions for P± in terms of the Laplace transform 
of Ay{t) and Ny{t) are readily obtained after the sub- 
stitution + ^ - and <<->>, and G(M) ^ G(M) in 
Eqs. (|0,|l|,|2l|6|J63). Finally, due to the symmetry of 



In terms of W{s) and V{s), the interval distributions take 
the symmetric form 



G{M) - sV{s) - sW(s) 
G{M)- sV{s) + sW{sy 



(76) 



p ^ G{M) + sV{s)-sW{s) 
G{M) + sV{s)+sW{s)' 
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In practice, the explicit forms of W{t) and V{t) obtained 
in Eqs. ( |67| , [70| ) permit a fast and efficient numerical im- 
plementation of the IIA. 

Note that for a smooth non- Gaussian process, all the 
above results of the IIA remain unaltered, G{M) now be- 
ing the cumulative sum of the associated distribution of 
X, and T being given by the general form po|. 



9{M){\X'{t)\), 



(78) 



whereas t± are still given by Eq. ([l5|). Applying the 
IIA results to a general non-Gaussian process only re- 
quires the prior knowledge of A< {t) and -/V< {t) (or A> {t) 
and Ny{t)). In general, these time-dependent functions 
should be given a priori, analytically, or extracted from 
numerical or experimental data. 

Finally, let us briefly address the validity of the IIA. 
First, crossing intervals are strictly never independent, 
except in the particular case of a Markovian process (see 
Eq. (^), for which the IIA does not apply, due to the 
singular nature of this process. The IIA is also an un- 
controlled approximation which seems almost impossible 
to improve systematically, by introducing interval corre- 
lations. However, in practice, the IIA is found numeri- 
cally to be a surprisingly good approximation, especially 
for "very smooth" processes for which f(t) is analytic 
|20|| (see the counterexample of /3(t) in Eq. (|90| ) be- 
low). We will even show in the next sections that some 
of the predictions of the IIA are in fact exact for smooth 
Gaussian processes. 



III. PERSISTENCE EXPONENTS 

A. General properties and physical applications 

In many contexts, one is interested in the large time 
behavior of the persistence probabilities P< (t) and P> (t) . 
It has been rigorously established that, if \f{t)\ < C/t, for 
sufficiently large time t (C is some arbitrary constant), 
then the persistence decays exponentially Hence, we 
define the two persistence exponents, by the asymptotics 



P<(t) 



P>{t) 



(79) 



valid when t — s- +oo. Due to the symmetry relation of 
Eq. (pa), we have 



(80) 



Hence, from now, we only consider the case M > 0. From 
Eqs. (16 17|Jlq , [l9| ) , we find that the interval distributions 
P± (t) decays in the same way as their associated persis- 
tence, for t — > +00. 

The name persistence "exponent" (instead of "decay 
rate") arises from its numerous applications in out of 
equilibrium physi cs joj |lO|, ^ ^ ^ 0, 

1^, H H, |2j||~|6f Indeed, in many cases, the 



normalized two-times correlation function of the relevant 
physical variable Z(T) obeys dynamical scaling, 



{Z{T)Z{T')) 



F{T/T'), 



(81) 



where T is the physical time, and F is the scaling corre- 
lation function. Defining 



i = lnr, X{t) 



Z{T) 



(82) 



the resulting process X{t) becomes stationary in the new 
effective time t 0| , with correlator 



/(0 = i^(expi). 



(83) 



Hence, the persistence Py(M,t) for the process X(t) is 
equal to the probability Py{M, T) that the process Z{T) 
remains above the level M x {Z^{T)) up to time T = 
exp t ijl^, |l^ . Since {M, t) decays exponentially, the 
persistence of the process Z{T) decays as a power- law, 
hence the name persistence "exponent" , 



P^{M,T) = P^{M,t)^e 
Similarly, one has 

P^{M,T) = P^{M,t) ~ e 



(84) 



(85) 



In particular, for M = 0, the persistence of the pro- 
cesses X and Z are both equal to the probability that 
the associated process does not change sign up to the 
time t — \nT . 

In order to illustrate the dynamical scaling of the cor- 
relator resulting from Eq. (^), let us give three physical 
examples. In the next sections, our analytical results will 
be tested on the corre lators /i [t] , f2it), and fait), intro- 
duced below in Eqs. 

• Consider a c?-dimensional ferromagnetic system (for 
instance, modelized by the Ising model) quenched 
from its equilibrium state above the critical tem- 
perature Tc, down to Tc (critical quench) or below 
Tc (subcritical quench). As time T proceeds, corre- 
lated domains of linear size L{T) ^ T-^/^ grow, and 
this coarsening dynamics leads to dynamical scaling 
for the total magnetization, of the form Eq. (|8l|). 
Initially the spins have only short range spatial cor- 
relations, and as the domains grow, the correlation 
length remains finite, of order L{T). If L{T) re- 
mains much smaller than the linear size of the sys- 
tem, the law of large numbers ensures that the mag- 
netization Z{T), which is the sum of the individual 
spins, is a true Gaussian variable. For M = 0, the 
persistence is equivalent to the probability that the 
magnetization never changes sign from the time of 
the quench (T = 0), up to time T. For a critical 
quench, the persistence decays as a power-law with 
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a persistence exponent 9c, which is a universal crit- 
ical exponent of spin systems, independent of the 
familiar ones (/3, 77, z,...), due to the non-Markovian 
nature of the magnetization |l6) . For a subcrit- 
ical quench, the magnetization persistence also de- 
cays as a power-law with a universal d-dependent 
persistence exponent controlled by a zero temper- 
ature fixed point ||l^ , and the dynamical exponent 
is z = 2. 

• If the field 2'(x, T) evolves according to the d- 
dimensional diffusion equation (or more sophisti- 
cated interface model equations ||l9|). 



dZ_ 
df 



(86) 



starting from an initial random configuration of 
zero mean, the process becomes Gaussian for large 
times, another consequence of the law of large num- 
bers. For any fixed x, the normalized two-times 
correlator of 2'(x,r) obeys dynamical scaling, and 
the probability that Z(x, T) does not change sign 
decays as a power-law with a d-dependent exponent 
computed approximately in ||l5| . The associated 
stationary correlator in the variable i = In T is 



hit) 



1 



cosh^ (I) 



(87) 



Moreover, in d = 1 and at a fixed time T, the 
process Z{x,T) is a stationary Gaussian process 
in the spatial variable x. The variable X{x) = 
Z{x, T)/y/{Z^{x,T))^ (where ( . )^ denotes the av- 
erage over the spatial variable x) has a Gaussian 
correlator. Hence, we shall later illustrate our re- 
sults using the correlator 



/2(t) = e--. 

• The random acceleration process 
its equation of motion. 



dT 



(88) 

is defined by 



(89) 



where t}{T) is a (5-correlated white noise. Again, its 
two-times correlator obeys dynamical scaling, and 
the associated stationary correlator is ll^ , 



1 3|t 

- e 2 
2 



(90) 



For M = 0, this process is a rare case for which the 
exact value of the persistence exponent is known 
exactly §, 



:(Af = 0) 



(91) 



Note that for the random acceleration process, the 
correlator f-i{t) is not analytic. Although, twice 



differentiable at t = 0, its third derivative is not 
defined at t = 0. For this process, it is not sur- 
prising to find that the IIA is not as precise as for 
smoother processes [|l^, Indeed, one finds the 
IIA resuh 6l±(M = 0) = 0.2647..., off by 6% com- 
pared to Eq. (^l|) , a relative error much bigger than 
usually observed for M = persistence exponents 
obtained by means of the IIA. 



B. Persistence exponents within the IIA 

Within the IIA, the persistence exponents 0± are ob- 
tained as the first pole on the negative real axis of the 
Laplace transform of the associated interval distribu- 
tion P±{t), since the Laplace transform of exp{—9±t) is 
l/(s + 6±). Using Eqs. (|6^,|6^, we find that 9± satisfies 



G(M)[1 + 0_iV<(-0_)] + 0_i<(-e_) 
G(M)[1 - 9+N^{-e+)] + 9+A^{-9+) 



1 



Tt 

0, 



or the equivalent relations in terms of and 
G{M)[l + 9+N>{-9+)]+9+A^{-9+) - ^ 
G(M)[1 - 6i_iV>(-6l_)] + 6i_i>(-6i_) 



r ( 
0. 



■,(92) 
(93) 

.,(94) 
(95) 



In terms of the auxiliary functions V{t) and W{t) intro- 
duced in Eqs. ( |6^ , |6^ , |6^ ) , the defining equations of 9± 
take a simpler form 



Wi-9+) - V{-9+) 



W{-9^) + V{-9^) = 



G{M) 



G{M) 



(96) 
(97) 



The residues R± associated to 9± , and defined by 

P±{t)^R±e-'^\ (98) 

can easily be extracted from Eqs. ( |7^ , |77| ), using the iden- 
tity 



dP^ 

ds 



(99) 



Before addressing the limit M ±00 in the next sec- 
tion, we present some numerical results for moderate M, 
and for the processes associated to the correlators fi{t), 
f2{t), and fsit), introduced in Eqs. ( ^7p^j90| ). We recaU 
the symmetry relation 9±{—M) = 9zp{M), so that we re- 
strain ourselves to the case M > 0. In order to compare 
the values of 9± for the different correlators, it is instruc- 
tive to multiply the persistence exponents (of dimension 
[t]^^) by the time scale tq, which is the mean crossing 
time interval for M = 0, 



To = t{M = 0) = 



(100) 
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We have performed extensive numerical simulations of 
the processes associated to the correlators fi{t), f2{t), 
and fsit), and measured the persistence and the crossing 
intervals distributions, and in particular, the persistence 
exponents. Let us briefly describe how to generate long 
trajectories of a stationary Gaussian process solely char- 
acterized by its two-times correlator p^ . In real time, 
the most general form of such a process X reads 



X{t) = f 



J{t~t')r]{t')dt', 



(101) 



where ri{t) is a (5-correlated Gaussian white noise. The 



Gaussian nature of 'r]{t) and the linear form of Eq. (101) 
ensures that X{t) is a Gaussian process. Moreover, the 
convolution product of the noise 'ri{t') with the kernel 
J{t — t') (instead of a general kernel J{t, t')) ensures sta- 
tionarity. Taking the Fourier transform of Eq. (101), we 
obtain 



X{u;)^J{u;)f,{Lo), 



(102) 



where X{ijj) — X{t)exp{—iujt) dt, and the noise 

Fourier transform satisfies {fj(uj) fjluj')) — 2t: 5{u} + lo'). 
The Fourier transform of the correlator of X is hence 



{X{u)X{co')) 
2n 



= f{uj) 6{UJ + Lo') = I J(c<j)|2 6{UJ + LJ'), 



(103) 

which relates the kernel J{t) to the correlation function 
/(i), through their Fourier transform, /(w) = |J(a;)p. 
Note that a necessary condition for f(t) to be the corre- 
lator of a Gaussian process is that its Fourier transform 
/(w) remains positive for all real frequencies lj. Finally, 
a trajectory of X is obtained after sampling 



(104) 



on a frequency mesh, and performing an inverse fast 
Fourier transform of the obtained X{bj). The Fourier 
transform of the correlators /i(<), /2(i), and f'i{t) having 
simple explicit expressions, this procedure for obtaining 
long trajectories is extremely efficient. Otherwise, one 
has to tabul ate t he Fourier transform of /(t), before sim- 
ulating Eq. ( p^ . 

In practice, our numerical results are obtained after 
averaging 10^-10^ trajectories of length T = 1024 tq 
(sometimes 2048 tq) and with a frequency mesh of typ- 
ically 1024^ points spaced by Atj = T^^. Since for a 
general M, r = roexp(M^/2), we obtain typically 10^~ 
10^*^ X exp(— M^/2) M-crossings, a number which decays 
rapidly for large M . Despite the loss of statistics for 
large positive M, we find that 9-{M) can still be mea- 
sured with great accuracy, since P- {t) and P< {t) becomes 
purely Poissonian in this limit (see next section). On the 
other hand, the error bars for 9+{M) increase rapidly 
with M, due to the occurrence of less M-crossings. In ad- 
dition, the determination of 6'+(M) is also plagued by the 
fact that the exponential asymptotics of P+ {t) and P> [t) 



M 




a sim 


T^O tl-,2 


n sim 





1.1700 


1.178(2) 


1.2928 


1.330(5) 


1/2 


0.6949 


0.7008(7) 


0.7587 


0.7723(8) 


1 


0.3715 


0.3743(6) 


0.3994 


0.4032(8) 


3/2 


0.1734 


0.1750(4) 


0.1831 


0.1850(4) 


2 


6.813 ■ 10-2 


6.834(6) ■ 10-2 


7.065 • 10-2 


7.116(7) ■ 10-2 


5/2 2.177 ■ 10-2 


2.180(3) ■ 10-2 


2.224 • 10-2 


2.225(3) ■ 10-2 


3 


5.509 ■ 10-^ 


5.510(2) ■ 10-^ 


5.568 • 10-3 


5.568(2) • 10-3 



TABLE I: We report the value of tq 6-{M) as obtained from 
the IIA calculation (OH:^) and simulations (61™), for different 
values of M, calculated for the processes associated to the 
correlators fi{t) (ro — 2n) and /2(f) (tq = tt), introduced in 
Eqs. (g,||). 



only develops for increasingly large t, as M increases, 
which may produce uncontrolled systematic errors in the 
numerical estimates of 6+{M) (see next section). 

Tabic I compares tq 6-{M) for the correlators fi{t) (for 
d = 2) and f2it), as obtained from numerical simulations 
and from the IIA calculation of Eq. (|9^) . The agreement 
between the theory and the simulations is excellent for 
both correlators, and is even improving as M increases. 
In fact, we will show in the next section that the IIA 
becomes exact for 9-{M), when M ^ 00, and we will 
obtain an analytic asymptotic expression for 9- (M) . Our 
theoretical and numerical results are also consistent with 
the numerical bounds computed in Q for the correlator 
/i(t), for M = 1 and AI = 2: 

0.3681 < To 6i_,i(Af = 1) < 0.4298, (105) 
0.0666 < To 6i_4 (M = 2) < 0.0748. (106) 

It is also clear that the IIA provides much better es- 
timates of the persistence exponent 9-{M) than these 
bounds, which are however exact, although they require 
a much bigger numerical effort than the IIA Q . 

Table H compares tq 9+ (M) for the correlators /i (t) 
and f2(t), as obtained from numerical simulations and 
from the IIA calculation of Eq. (96). The agreement 
between the theory and the simulations is satisfying, al- 
though we will show that the IIA ultimately fails in pre- 
dicting the exact asymptotics of 9^(M) as M ^ 00, al- 
though the limiting form of P+{t) will be given exactly 
by the IIA, for t not too large. 

For the process associated to the correlator /a (t) (tq — 
2tt/V3), to9^^^{M = 0) = 0.9602..., compared to the 

0.9069... obtained 



0) 



exact value tq 9± {M — — 

in [0. The agreement between the IIA and numerical 
estimates of 9^ [M) greatly improves as M increases, as 
observed for the two other correlators in Table |, whereas 
9+{M) is only fairly reproduced for large M (as already 
observed in Table y) , suggesting that the IIA somewhat 
fails to reproduces 9+{M) in the limit M -l-oo, as 
will be confirmed in the next section. Finally, we note 
that for the three processes considered here (and all other 
smooth Gaussian processes known to the author) , one has 
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^ /iIIA ^ /I sim 



1.1700 1.178(2) 
1/2 1.8164 1.865(6) 

1 2.6475 2.67(2) 
3/2 3.6677 3.74(3) 

2 4.8926 5.06(5) 
5/2 6.2651 6.6(1) 



1.2928 1.330(5) 
2.0232 2.125(7) 
2.9606 3.11(2) 
4.1031 4.36(4) 
5.4338 5.90(7) 
6.9152 7.6(2) 



TABLE II: We report the value of ro 6+{M) as obtained from 
the IIA calculation (S+*) and simulations (0+™), for different 
values of M, calculated for the processes associated to the 
correlators fi{t) (ro = 2-n) and /2(f) (to = n), introduced in 
Eqs. (pTia. 



B±{M = 0) - Tq^S with a proportionality constant close 
to unity. 



Exact results in the limit M — > ±00 



In the limit M — s- +00, and using 

_ P 2 

G{M) 



/27r M 

and Eqs. (|l^,|l^), the mean of the ± intervals are 

2lT A: 



: e 2 



2r, 



(107) 



(108) 



(109) 



Hence, as expected physically, the typical length of the — 
intervals is becoming increasingly large as M — > +00 (of 
order 2t), whereas the typical length of the + is going 
slowly to zero. 



1. Distribution of — intervals and 6- 

For large M, and using Eqs. ( |67|j70| ), we obtain that 
V{t) ^W{t)(l + o(e-^)) . (110) 



Moreover, in this limit both functions can be approxi- 
mated by developing f{t) up to second order in t, 



1 M^ant 

V{t)^W{t)^-e 



1 



(111) 



The constraint that 

^+00 



.-1 







W{t)dt = A<(0) - A<(+oo), (112) 

= A>(0)-A>(+cx3), (113) 
= G{M)G{M), (114) 

a 9 

(115) 



'2tt M 



is consistently recov ered up to leading order, after in- 
tegrating Eq. (111). Moreover, the explicit form of 
Eq. (Ill) implies that the Laplace transform of W{t) 
can be approximated by its value at s = 0, provided that 



M> 1. 



(116) 



Finally, using Eqs. (^,^^, we find that the Laplace 
transform of P_ [t) is given by 



1 



1 



(117) 



under the conditions of Eq. ( |116| ). Hence, we conclude 
that for large M, one has 



(118) 



and, taking the inverse Laplace transform of Eq. (117) 
that the distribution of — intervals is essentially Poisso 
nian 



P-{t)=TZ^e 



P<{t) 



(119) 



except in a narrow region of time < i < T4. <C t_, 
correspondin g to the conjugate time domain of the con- 
dition of Eq. (116). The behavior of P-{t) and P<{t), for 
t ^ T+ , will be obtained exactly in section [V . The fact 
that the distribution of the long — intervals is becoming 
Poissonian can be physically interpreted. Indeed, since 
T_ +00 when M +00, the process X can be con- 
sidered to be Markovian at this time scale. In addition, 
this also shows that for the — intervals, the IIA is in fact 
exact. In the next section, in the process of finding exact 
bounds for 9+{M), we will prove the exact result, valid 
in the opposite M ^ limit. 



2/(0) 



M + O {M^ 



(120) 



where 6*0 = 0_{M = 0), {X{t))x>o is the average of X 
over all trajectories for which X{t) > for all times, and 

On Fig. 0, we plot 6-{M) obtained by simulating the 
process associated to the correlator /2(i) of Eq. (|8^). We 
also plot the IIA result of Eq. (p7|), which is in perfect 
agreement with numerical simulations. In addition, we 
illustrate the rapid conv ergence of 0_(M) to the exact 
asymptotics of Eq. (118). On Fig. |[ we plot P-{t) for 
M = 3, which has already perfectly converged to its 
asymptotic Poissonian form of Eq. (lllJ 



2. Distribution of 4- intervals and 0+ 

For M +00, and using the asymptotics Eq. ( |110| ) 
and the IIA expression for P+{s) of Eq. (fz^), we obtain. 



P+{s) = 1 -2sW{s). 



(121) 
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FIG. 1: (Color online) We plot 6-{M) (in unit of r,^^), for 
the non-Markovian process associated to the correlator f2{t) 
of Eq. (j^j . Symbo ls ar e the results of numerical simulations 
(see text and Eq. (104)), while the full line is the result of 
the IIA approximation, Eq. (p7|). Finally, the dashed line 
corresponds to the exact asymptotic result, d-{M) ^ tZ^ ■ 
Some values of 6- (Af ) are also reported on Table ^ 




FIG. 

of e 



2: (Color online) For M = 3, we plot the distribution 
— ± intervals, P±{t), as a function of t/T± (full lines; 



for M = 3, T+ « 0.76350 a. 



-1/2 



and 



564.83 a. 



-1/2. 



for 

the process associated to the correlator /2 (t) of Eq. (M) . The 
distribution of — intervals has converged to the Poissonian 
form of Eq. (119) (dashed line). N ote that the linear regime 



168), for f <C r+ ^ r_, cannot 



of P- {t) predicted in Eqs. ( p^ , 
be seen at this scale. The distribution of + inte rvals is well de- 
scribed by the Wigner distribution of Eq. (122) (dashed line), 

(insert; 



but ultimately decays exponentially for t 3> 
dotted line of slope r+ 9+ ^ 2.35). 



-1/2 



In real time, Eq. (121) reads 

p+{t) ^ -2TA';{t) ^ 



TT t -1 

2rl 



(122) 



which is vahd for t <ti . Using Eqs. (23 10§| ), the 



relation P+{t) = — 2tA" (t) can be explicitly written as 

P+it) = -T^{X'{t)X'{0) S{X{t) - M) 6{X{{)) - M)). 

(123) 



When the process X crosses the level M for the first time 
at time after the preceding crossing at time 0, one has 
-X'{t)X'{fi) > 0. Eq. (|T2^) states that for small times 



— 1/2 

t <C , the probability of having a + interval of length 
t is of the same order as the probability of having an M- 
crossing at time t, knowing that there is such a crossing 
at time 0, which follows a — interval (hence the factor 
T_). In other words, for short time a crossing at time 
t is very often the first crossing following the crossing at 
time 0. 



The form of P+{t) found in Eq. (122) can be obtained 
from a simple physical argument, which will be our basis 
for obtaining exact results for the small t behavior of 
P±{t) in section IV. Indeed, let us consider a + interval 

]^ /2 

of length < ^ flj . For < t' < t, one can expand the 
smooth process X{t') in power of t' up to second order, 
starting from X{0) = M, 



X{t') 



M - 



X'iO) t' + ^^t'' + 



(124) 



remembering that X'{t) is a Gaussian variable of second 
moment (X'^(t)) = 02 = -/"(O), the probability distri- 
bution of the velocity v = X'{0) > 0, at a M-crossing 
following a — interval, is given by 



p{v) 



(|X'(0)|<5(X(0) - M) SiX'iO) - z.))x^(o)>o 



(|X'(0)|5(X(0)-M))^,(„)>o 

(125) 

Since X{0) and -'^'(O) are uncorrelated (as 
(X'(O)X(O)) = /'(O) = 0), these Gaussian averages 
are easily performed, leading to the exact result 



p{v) = —e~ 

0.2 



(126) 



In addition, the distribution of X"{0), conditioned to the 
fact that X{0) = M, is a Gaussian of mean 



{X"iO))xio)=M = -Ma2, 



(127) 



and mean square deviation f^^{0) — /"^(O), which is of 
order a^, and is independent of M. Hence, for large M, 
one can replace X"{0) by its average, and the interval 
length t can be obtained by finding the first M-crossing 
of the trajectory of Eq. (124), 



2v 



Ma2 

Hence, for small time t, on has 



P+it) 



+00 



p{v) 5{t — ty) dv. 



(128) 



(129) 



Finally, using Eq. ( 126 ) and the asymptotic expression 
for T+ of Eq. ( |108| ), we o btain the distribution of + in- 
tervals given by Eq. ( |l22|) . 

Thus we find that the distribution of + intervals is 

— 1/2 

given by the Wigner distribution for t , at re- 

sult also obtained in [01. However, note that the ratio 



12 



t/rj. can be arbitrarily large in the limit M — > +cx3, or 
T_|_ 0. The probability distribution of Eq. (122) is cor- 



rectly normalized to unity and has a mean equal to t+. 

— 1/2 

Of course, for t , the actual distribution of + in- 

terval should decay exponentially 1^, with a rate 9+{M). 

— 1/2 

Matching the two asymptotics at t ~ 02 , we find that 
up to a so far unknown multiplicative constant, 



(130) 



for large M . In addition, the above argument shows that 
the total probability contained in the exponential tail of 
P+{t) vanishes extremely rapidly as M ^ -|-oo, and is 
of order exp(— iiTM^), where K i s a c onstant of order 
unity. Moreover, the result of Eq. (122) implies that the 
persistence is given by 



. 4 ( x+ ) 



1/2 . 1/2 

for t <C 02 , and decays exponentially for t ^ a2 

In the limit M +00, the determination of d^{M) 

seems to be beyond the IIA. Indeed, let us assume for 

simplicity that f(t) decays faster than any exponential, 

f{t) ^ exp(— en), with 7 > 1 . Anticipating that 9^{M) 

is large, we need to evaluate P+{s) for large negative s. 

In this limit, Eqs. (0j7|) leads to 



V{s) 
W{s) 



M 



2tt 



e 2 s/(s). 



(132) 
(133) 



If 1 < 7 < 2, and using the IIA expression for P+{s) of 
Eq. ([76|), we obtain 



(134) 



which leads to 



5IIA 



(M) 



2 02 



(135) 



which grossly underestimate the divergence of 9^{M), 
when AI — > 00. If f{t) ^ exp(-'Ct^) decays as a Gaussian 
(7 = 2), Eq. (M) leads to 



qIIA 



(M) 



(136) 



which, again, behaves linearly with M. Finally, if 7 > 2, 
we find that 



(137) 



up to a computable multiplicative constant. 

Let us now present exact bounds for P> {t) which will 
lead to an exact asymptotics for 6'_|_(Af), fully consistent 



with Eq. ( |130| ). We discretize time U = iAt, with At 
t/n, and define Xi — X(ti). By definition. 



/•+00 

P>{t,M) = 

J M 



Im (27^)"/2^/detC 
where the Gaussian "action" has the quadratic form 



(138) 



(139) 



and where the matrix D is the inverse of the correlation 
matrix C defined by its matrix elements 



{X{U)X{t,))=f{U-t,). 



(140) 



Making the change of variables Ui — Xi + M E [0, +00], 
and noting that 



(131) si{y,}) = Iy1 y^yj + E y^ + ^H ' 



2J 



we obtain 
P> {t, M) = P> [t, M = 0) {^e^^'^S. 

with 



n 1 " 



(141) 



(142) 



(143) 



and where ( . ) y>o denotes the average over all processes 
for which y{ti) > 0, for all i. If we assume that the 
process is periodic of period t, the vector u = (1, 1, 1) 
is an exact eigenvector of C associated to the eigenvalue 



n/2 

A= E /(^^)' 

i=-n/2 



and one has 



(144) 



(145) 



For large time t, the periodic constraint should not af- 
fect the value of the persistence exponent. In fact, in 
this limit of large time and fine discretization {At 0), 
and even dropping the assumption that the process is 
periodic, one finds that 



At 



At 



jr^mdt 2/(0)' 



(146) 



where the discrete sum of Eq. (144) has been transformed 
into an integral when At 0, and the integral limits 
±t/2 extended to icxi for large t. Finally, defining Oq = 
6± (M — 0) and using Eq. ( 142 ) , we find the exact result 



0+{M) = 60 + e{M) + 



47(0)^ 



(147) 
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with 



e(A/) 



lim — In ( e ^ 

t^+oo t 



■M- r xit')dt' 

/(O) Jo ^ ' 



(148) 



X>0 



Using the convexity of the exponential function, and the 
fact that the argument in the exponential in Eq. (148) is 
negative, we obtain the exact bounds 



0<e(M)<^^^2^M, 
" 2/(0) 



(149) 



where {X{t))x>o is the average of the process X re- 
stricted to the trajectories which remain positive for 
all ti mes, and is a constant strictly independent of M. 
Eqs. (|147| , |149|) lead to an exact bound for 0+{M), 



eo + ^<e+ (M) <9o + i^^J^^ M - 

4/(0) - 2/(0) 
which implies that for large M, 

9+{M) 



M 



,t2 



M2 



4/(0) 



4/(0) 
(150) 



(151) 



with a subleading correction bounded by a linear term 
in M . The exact asymptot ics o f Eq. (151) confirms our 
heuristic argument of E q. ([130^. In addition, for small 
M, it is clear from Eqs. ( |148| , p-49( ) that one has the exact 
expansion. 



{X{t)) 



M + O (M^) . 



2/(0) 

Since 6i_(M) = 61+ (-M), we also get 

2/(0) ^ ' 



(152) 



(153) 



On Fig. H, we plot 6'+(M) obtained by simulating the 
process associated to the correlator f2{t) of Eq. (88). 
We also plot the IIA result of Eq. (^), which under- 
estimates the actual value of 0+(M), as explaine d ab ove. 
In addition, we plot the exact bounds of Eq. ( |150| ), as 
well as a convincing fit of 9+ (M) , to the functional form 
9+{M) = ao + ai M where the exact leading term 



was obtained in Eq. (151). On Fig. ||, we plot P+{t) for 

M = 3, which follows the predicted Wigner distribution 

n K —1/2 

of Eq. (122), for t ^ flj , before decaying exponen- 
tially, with rate 6'+(M), for large t. 




FIG. 3: (Color online) We plot e+{M) (in unit of t^^), 
for the non-Markovian process associated to the correlator 
f2(t) of Eq. (^8|). Symbols correspond to results of numerical 
simulations. T he u pper and lower dashed line are the exact 
bounds of Eq. (150) (the upper bound being exact up to order 
M, for small M). The ratio of these two bounds goes to 1 
when M — > +CX). The full line is a quadratic fit for M > 3/2, 
to the functional form 6j^{M) = ao -|- Oi M -|- jj^j with oq « 
2.0 and ai ~ 0.66. Finally, the middle dashed line is th e IIA 
result (along with its asymptotic slope given by Eq. ( |l36[ ); 
dotted line), which underestimates the quadratic growth of 
9+{M). Some values of 0+{M) are also reported on Table ^ 



differentiable, we will show that the IIA surprisingly re- 
produces these exact results. However, for a marginally 
smooth process, such that for small i, and 2 < a < 4, 



(154) 



the fourth derivative of fit) is not defined at 0. We will 
show that for such a process, the IIA does not lead to ex- 
act results for the small t behavior of P±{t), although it 
is in fact qualitatively and even quantitatively accurate. 
Note that the process associated to the correlator fs(t ) 
defined in Eq. (90) satisfies the property of Eq. ( 154 ), 
with a = 3, which implies that its velocity is not differ- 
entiable. 



A. Exact small time behavior for very smooth 
processes 

In the limit t ^ 0, the trajectory of the very smooth 
process X inside a -I- or — interval is essentially parabolic. 



IV. DISTRIBUTIONS P±(t) FOR SMALL 
INTERVALS 



The heuristic argument presented in the preceding sec- 
tion, below Eq. (124), can be adapted to provide the ex- 
act behavior of P± (t) , P>{t), and P<{t) for small time 
t. For a smooth process with f{t) at least four times 



X{t') = M + vt' + ^t'^ + 0{t^), 



(155) 



where the distribution of the velocity at a crossing time 
is given by 



p{v) = —e~ 
a2 



(156) 
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and the distribution of the acceleration a, conditional to 
the fact that X{0) = M, is 



cr(a) 



1 



27r 022: 



with 



/"2(0) 



-1, 



(157) 



(158) 



and a2 — ^/"(O)- Note that the acceleration at 
t' = is independent of the velocity at t' = 0, since 
{X"{0)X'{0)) = /'"(O) = 0. In addition, since 

/^^(O) - r (0) = {u' a,y m ^ > 0, 

(159) 

z defined by Eq. (158) is indeed a positive real number. 



Let us first consider the small time behavior of P+ (t) 
For an interval of length t to be small, and since v > 
the acceleration is necessarily negative. From Eq. (155) 
the crossing time is then given by 



t ■ 



2v 

a 



(160) 



which is only valid when t is small. Hence, for small t, 



P+it) 



+00 

dv 

J -00 



da piv)<7{a)S [t + — ] . (161) 



After integrating over u, we obtain 



(162) 



Using the explicit form of p(w), and taking the limit t —> 
0, we obtain 

P+{t)=c+{M)a2t + 0{t^) , (163) 
where the dimensionless constant c+(M) is given by 

1 

c^(M) = a{a)a^da. (164) 

4 0'2 J -00 

Performing explicitly the Gaussian integral above, we fi- 
nally obtain the exact result. 



zM f M\ , , 
— .9 -), (165) 



where g{X) is the Gaussian distribution and GiX) its 
cumulative sum, both defined in Eqs. (^|J^). In the limit 
M — > -l-cxD, we find c+(A'f) - M2/4, which leads to 



(166) 



Oh 



"0.5 




FIG. 4: (Color online) For M = 1/2, we plot the distri- 
bution of £ = ± intervals, P±{t), as a function of t/T± (full 
lines ), f or the process associated to the correlator f2{t) of 
Eq. (BSl) (P- (t) is the most peaked distribution) . The straight 
dotte d lines have the predicted slopes at f = 0, given by 
Eqs. (163 165, The dashed lines are the distributions 
rom the IIA, after taking the inverse Laplace trans- 
(|7q,[r7|). The insert shows the same data on a 



obtained 
form of Eqs. 

semi-log scale, illustrating the good accuracy of the IIA in 
predicting the persistence exponents 6± (M) and their associ- 
ated residue R±. 



in agreement with our result of Eq. (122). Moreover, 
using Eq. (16), we obtain the small time expansion of 
P>(t), up to third order in time, 

P->it)^l-^ + ^t' + 0{t'). (167) 



6t+ 



Finally, the corresponding results for c_(Af), P-{t), and 
P<{t) are obtained by the substitution M ^ ~M, and 
the exchange of the indices -!-<-)■ — and > ^ < . In par- 
ticular, we have 



c_(M) 



A/2 -I- ^2 



G 



M 



z 



zM f M 

—r9 



In the limit M 



00, we find that 

z^ /M 
2M3^ I T 



(168) 



(169) 



Thus, we find that P±(t) behaves linearly with time for 
small t, a result which will be shown in section IV C to be 



specific to very smooth processes, for which the correlator 
f{t) is at least four times differentiable at t = 0. 

On Fig. |, for M = 1/2, we plot P±{t) obtained 
from numerical simulations of the process associated 
to the correlator f2{t), illustrating their linear behav- 
ior for small time t. The initial slope at i = is in 
perfe ct numerical agreement with the exact results of 
Eqs. (163 165 |l6^ ). In addition, we also plot the full dis- 
tributions P± {t) , obtained by taking the inverse Laplace 
transform of Eqs. (|76| , [77| ). For this moderate value of 
M, far from the large M regime where the IIA becomes 
exact, the good agreement between the IIA results and 
numerical simulations is certainly encouraging. 
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B. II A results 



fait) introduced in Eq. (Q), we find that P±{t) should 
be constant at t — 0, with 



We now derive the small time behavior of P±{t) by 
using the IIA. Expanding the explicit form of V{t) an 
Wit) of Eqs. we find, 



as 
0-2' 



(177) 



1 M2 4- ^2 

Wit) = i_i!i-±^a2t2 + 0(i4^ 



Vit) = GiM) 



16 
1 



(170) 



G 



M 



For the correlator fsit), one has 02 = 3/4 and 03 = 1/4, 
so that P±(0) = 1/3. 

Let us now derive an exact expression for P+it) for a 
marginally smooth process with a = 3. For small t, the 
correlator of the velocity is 



zM 



Taking the Laplace transform of these expression and in- 
serting them into the IIA expression for P+is) of Eq. (76) 
we obtain, for large s. 



8 \ \ z 

a^t^ + Oit^). (171) {X\t)X\Q))^-fit) = a2-6a3\t\ + 0{t^), (178) 

which coincides with the small time behavior of the cor- 
relator of a Markovian process (see Eq. (0)). Hence, on 
short time periods, the local equation of motion of X'it) 
is. 



P+is) 



+ zM g 



(M^ + z^) G 



M 



M 



0-2 

4s2 



+ o(.-4). 



(172) 



which is exactly the Laplace transform of Eqs. (163,165). 
Thus, we find that the IIA reproduces the exact small 
time behavior of P± it) . 



C. Marginally smooth processes 

In this section, we study marginally smooth processes 
characterized by a correlator fit) having a small time 
expansion of the form of Eq. (154), so that /^^(O) does 
not exist. Since /^^(O) was explicitly appearing in our 
results of the preceding sections, this suggests that the 
small time behavior of P±it) should be affected by the 
weak singularity in fit). 

Let us first apply the IIA in this marginal case. The 
small time expansion of Vit) and Wit) now read, 



Wit) ^ i-^^^r- 

2 2 02 



(173) 



Vit) ^ G(M)- (174) 

Note that since 2 < a < 4, one has a — 2 < a/2. For 
large s, and using Eq. ( [76| ) and the above asymptotics, 
we obtain. 



P±is) 



(a — l)r(Q: — l)aa 



2a2s"-2 

In real time, this leads to the small time behavior, 
(a - l)r(a- l)aa 3 



P±it) 



2r(a-2) 02 



(175) 



(176) 



X"it) 



603 

0,2 



X'it) + 2V3^ vit), 



(179) 



where rjit) is a J-correlate d Ga ussian noise. From 
the equa tion of motion Eq. (179), one indeed recovers 
Eq. (ITS), for short time. Now using Eq. (179), a short 
time trajectory of X, starting from X(0) = M and 
X'iO) = V, takes the form. 



Xit) = M + V t + 2^/3^ Zit) + O it^ 



(180) 



where Z(t) — O (^'^'^ is the random acceleration process 
introduced in Eq. (|89|), 



Zit) = f dti [ dt2Vit2). 

Jo Jo 



(181) 



Finally, for small t, the first M-crossing of the process 
X, corresponds to the first time for which Zit)/t crosses 
the level 



Zo 



2x/3^ 



(182) 



Introducing the probability distribution "^ito.Zo) that 
Zit)/t crosses Zq for the first time at time t — to, one 
has the scaling relation 



^ito,Zo) = -2^ 



72 



(183) 



obtained by noticing that the scale invariant process Zit) 
has dimension [i]'^/^. For small time t, 



r+00 

P±it)= , 
Jo 



12a3t\ 12 as 



■ dv. 



(184) 



After making the change of variable T = 12 a^t/v^ and 
taking the limit i — > 0, while using the fact that p(w) ~ 
v/a2 for small v, we obtain the final exact result. 



which is independent of M . In particular, for the quite 
common case a = 3, which corresponds to the correlator 
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Up to a dimensional constant 603/02 depending on f(t), 
P-l-(O) is proportional to the mean inverse first-passage 
time of the process Z{t)/t at the level Zq = 1. By simu- 
lating the process Z, we have obtained 



Finally, within this simple approximation, we find that 

(l/T) = 1/3, overestimating the value obtained in 

Eq. PI). 



-)«0.193(1), 



which leads to 



P±(0) = 1.158(6) 



03 
02' 



(186) 



(187) 



This result has also been checked numerically for the pro- 
cess as socia ted to f3{t). In fact, the constant appearing 
in Eq. ( |187| ) was obtained exactly by Wong ||] , based on 
the study of the process Z{t) of Their result leads to 

P±(0) = g^ = 1.15625 ^. (188) 
32 a2 02 



Hence, we find that the IIA result of Eq. (|177| ) is not ex- 
act for marginally smooth processes with a = 3, although 
it predicts correctly that P± (0) is a constant independent 
of AI, leading to a reasonably accurate estimate of this 
constant. For general a, Eq. (176) is certainly correct 
dimensionally speaking, and probably fairly accurate, in 
practice. 

We end this section by an approximate calculation of 
{l/T) for the process Z{t), which does not reproduce the 
exact result of j^, obtained by a much more complex 
method O, m . We make the approximation. 



Zit) 



dU 



dt2'n{t2) 



t 
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dtniiti), (189) 



where the factor 1 / -\/3 ensures that both processes have 
the same mean square displacement {Z'^{t)) = t^/?,. 
Then, the original first-passage problem for Z{t)/t be- 
comes a standard first-passage problem for the usual 
Brownian motion B{t) at the level Bq = v^, for which 
the first-passage time probability distribution is given by 
ill, 



V. CONCLUSION 

In this work, we have considered the M-crossing in- 
terval distributions and the persistence of a smooth non- 
Markovian Gaussian process. We have obtained exact 
results for the persistence exponents in the limit of a 
large crossing level M, including exact bounds for 9+{M) 
and P^{t). In this limit, we have shown that the distri- 
butions of -|- and — intervals become universal, and are 
respectively given by the Wigner and Poisson distribu- 
tions. For any value of A/, we have obtained the exact 
small time behavior of the interval distributions and the 
persistence. We have also derived these results within 
the Independent Interval Approximation. Quite surpris- 
ingly, the IIA reproduces all these exact results, except 
for the large M asymptotics of 9+(M). In addition, the 
IIA fails in reproducing the exact small time behavior 
of the interval distributions for marginally smooth pro- 
cesses, although it remains qualitatively correct and even 
quite accurate in this case. To the credit of the IIA, it is 
the only method to provide precise approximate expres- 
sions of the interval distributions and the persistence for 
all times, and all values of the level M, and thus to grant 
access to the distribution of extrema of a non-Markovian 
Gaussian process. In addition, the IIA can be straight- 
forwardly applied to any smooth non- Gaussian process, 
for which the autocorrelation function A>(t) (or A<(t)) 
and the conditional number of crossings A^> (t) (or A^< (t)) 
are known analytically, or extracted from experimental 
or numerical data. For Gaussian processes, simple forms 
of the derivative of these quantities have been obtained, 
which permits a simple and fast numerical implementa- 
tion of the IIA results. 



Bn 



and for which, 



/2^T3/2 



B2 ■ 

^0 



(190) 



(191) 
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